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Abstract. Recent experimental work has suggested that electroosmotic flows (EOF) 
through conical nanopores exhibit rectification in the opposite sense to the well-studied 
effect of ionic current rectification. A positive bias voltage generates large EOF and 
small current, while negative voltages generate small EOF and large current. Here 
we systematically investigate this effect using finite-element simulations. We find 
that inside the pore, the electric field and salt concentration are inversely correlated, 
which leads to the inverse relationship between the magnitudes of EOF and current. 
Rectification occurs when the pore is driven into states characterized by different salt 
concentrations depending on the sign of the voltage. The mechanism responsible for 
this behaviour is concentration polarization, which requires the pore to exhibit the 
properties of permselectivity and asymmetry. 
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1. Introduction 

Nanopores are apertures with dimensions of a few tens to a few hundreds of nanometres 
in an insulating matrix. Immersed in an electrolyte solution, they can conduct a flux 
of ions and water through their interior. Modulations of this flux by electrically-driven 
translocation of charged macromolecules form the basis of single-molecule resistive pulse 
sensing, a major technological application of nanopores today DEI- 

Nanopores and nanochannels are characterized by large surface-area-to-volume 
ratios, and thus the transport of water and ions through these confined geometries 
is strongly dependent on their surface properties In particular, nanopores made 

from silica or glass typically develop a negative surface charge at neutral pH [5]. When 
in contact with a salt solution, the surface is screened by an electric double layer a few 
nanometres thick containing mobile positive counterions. Application of a tangential 
electric field results in an electric force on this screening layer. Viscous transfer of 
momentum away from the surface leads to bulk fluid motion known as electroosmotic 
flow (EOF) |6j. 

In addition to fluid flow, the screening layer affects ionic transport significantly 
as the dimensions of the system are reduced. As the conductivity due to the double 
layer becomes comparable to the conductivity of the pore’s lumen, the pore exhibits 
permselectivity towards cations HI IB]- 

In this regime of surface-governed transport, interesting effects have been observed. 
As with permselective ion-exchange membranes, driving an ionic current through a 
permselective pore results in concentration polarization (CP) across the pore [9]. In 
order for flux balance to be achieved, the total salt concentration reduces upstream 
of the permselective region, and increases downstream HBj. As the driving voltage is 
increased, the upstream depletion becomes more pronounced until the supply of charge 
carriers is exhausted, leading to a limiting current behaviour. At higher voltages still, 
clectrokinetic instabilities result in vortices which destroy the depletion region, and 
conduction is recovered through the pore mm- 
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If an asymmetry is introduced along the axial direction of the pore, for instance by 
making the pore conical in structure, the conductive properties of the pore now change 
with respect to voltage reversal. This is because one polarity of the applied voltage 
results in a depletion region inside the pore, while the other results in enrichment. 
The pore therefore exhibits ionic current rectification, and behaves like an ionic diode 

mmmm- 

The ionic transport properties of rectifying nanopores have been well-studied. 
However, in addition to ionic current rectification, there has been recent experimental 
evidence suggesting that electroosmotic flow (EOF) through conical nanopores also 
exhibits interesting and partly unexpected behaviour [T6J [T7, T8] . In particular, EOF 
has been observed to exhibit flow rectification which behaves in the opposite sense 
to current rectification, i.e. large currents are correlated with small EOF, and vice 
versa. The relative scarcity of studies on EOF rectification in single nanopore systems 
is most likely due to difficulties in measuring the effect, as compared with ionic current 
measurements. 

Rectification of EOF has potentially significant consequences as it could guide the 
development of nanofluidic flow rectifiers. These would play the same role in microfluidic 
circuits as electronic diodes do in electrical circuits. Due to the linearity, and hence 
time-reversal symmetry, of the Stokes equations governing low Reynolds number fluid 
flow, engineering low Reynolds number flow rectifiers has traditionally proven to be a 
challenging problem [12, ZD]- 

In this paper, we present the results of Unite-element simulations which capture 
the rectification behaviour of electroosmotic flows in a nanopore. We show that the 
fluxes of water and ionic current through the pore depend on the local electric held 
and salt concentration. By comparing our results with simple analytic models, we 
can predict the relationship between held and concentration inside the pore. Finally, 
we demonstrate how concentration polarization leads to variations in these quantities, 
resulting in rectification of EOF and ionic current which are intricately linked. 
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2. Model and Methods 

Finite-element simulations have proven to be a powerful technique in the analysis 
of nanopore physics [HI 121) 221 221 • For pores with the smallest dimensions larger 
than a few nm, the physics of ion and fluid transport is well-captured by classical 
continuum equations. Compared to atomstic and molecular dynamics simulations, 
continuum simulations can deliver accurate, quantitative results at a fraction of the 
computational cost, as long as the pore dimensions are sufficiently large that coarse- 
graining of molecular details is possible. We use the commercial software package 
COMSOL Multiphysics 4.4 to carry out finite-element calculations. 

Our simulation geometry is shown in Figure la. It is a 2D-axisynmretric geometry 
which can be revolved about the axis AH to produce full 3D solutions. The pore has a 
radius a = 7.5 nm and is embedded in a membrane of length L = 100 nm. The taper of 
the pore is parameterized by an angle a. The simulation is enclosed within a box of size 
d ~ 10L, which is large enough to mimic infinity (Supplementary Information). The 
membrane is modelled as a solid with permittivity e r ~ 4.2, similar to that of quartz 
glass; the solution is water with permittivity e r ~ 80. The surface of the membrane 
contains a constant surface charge density of a s . The constant a s approximation is 
reasonable at moderate salt concentrations (~ 100 mM) |24j . 

We use a free triangular mesh with a nine-level boundary layer at the surface 
CDEF. The thickness of the first boundary layer is set to ~ Ad/10, where Ad is the 
Debye screening length, to ensure the electric double layer is adequately resolved. The 
Debye length is the characteristic thickness of the electric double layer [3]. The mesh is 
heavily refined inside the pore as well as near the corners C, D, E, and F (Figure la). 

Transport through nanopores is well-described within the Poisson-Nernst-Planck- 
Stokes (PNPS) formalism [3]. The electric potential </>(r) is related to the charge density 
p e (r) by the Poisson equation: 

Pe(r) 


V 2 0(r) 


( 1 ) 
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Figure 1 . Finite-element simulations of transport through a nanopore, (a) The 
simulation geometry is a 2D axisymmetric design about the axis AH. The pore is 
modelled as a cone of angle a and tip diameter 2 a embedded in a membrane of 
thickness L whose surface CDEF contains a charge density a s . The box size d is 
chosen such that it mimics infinitely large reservoirs. A voltage is applied between 
two electrodes AB and HG. (b) Typical flow profiles through a pore with a = 0.1 
rad for an applied voltage of Vo = +1 V. Arrows indicate flow direction, and colours 
indicate flow magnitude, (c) Typical ionic flux through the pore for the same voltage. 
Arrows indicate the direction of net ionic current, while colours represent the total salt 
concentration. The bulk salt concentration is 2 cq = 200 mM. 























6 


where 6 q ~ 8.85 x 10 12 Fm 1 is the permittivity of free space. The ionic flux Ji of 


species i is calculated using the Nernst-Planck equation: 



( 2 ) 


where Di, Cj(r), and z t are the diffusion constant, concentration, and valency of species 
i respectively, e ~ 1.6 x 10~ 19 C is the elementary charge, ks ~ 1.38 x 10 -23 JK” 1 is the 
Boltzmann constant, T = 300 K the temperature, and u(r) the fluid velocity held. For 
KC1, Zi = ±1 and D K + ~ Dei- ~2x 10 -9 m 2 s _1 . Finally, the fluid velocity is obtained 
using the Stokes equation: 


fi'V 2 u = Vp + p e (r)V0(r), 


(3) 


where p is the dynamic viscosity and p the pressure inside the fluid. 

Our simulation proceeds in two steps. First, a ID solution is obtained along the 
boundary BCFG, with the boundary conditions (ft = 0 at B, (ft = Vo at G, surface charge 
cr s at C and F, and c t = Co at B and G. This results in solutions for the electric potential 
(ftiD and concentrations c i: in for both species along the boundary. The pressure p\r> can 
be obtained from the Stokes equation by substituting in the solution for (ftm. 

In the second step, the PNPS equations are solved fully coupled in the 2D geometry 
using the boundary conditions (ft = 0 on AB, (ft = V 0 on HG, surface charge a s on CDEF, 
and on AB and HG, no fluid stress, and q = Co- The solutions <ft\ d, Ci^d, and p±D are 
used as Dirichlet boundary conditions on BCFG. Our method is very similar to that 
used in a previous study [2TJ to model induced-charge electroosmosis through nanopores. 

We run the simulations for cq = 100 mM, and varying Vo, cr s , and a. The solutions 
are a set of fields <ft{r ), Cj(r), u(r), and p{r). Typical flow fields and concentration 
profiles are shown in Figure lb and c. In order to obtain the ionic current / and flow 
rate Q , we carry out the following integrals: 



(4) 



u ■ dS, 


(5) 
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where S is the surface area spanning the cross-section of the pore. For our calcuations, 
we chose the location of S to be halfway through the pore. Full details of the simulation 
procedure are given in the Supplementary Information. 


3. Results 


3.1. Rectification of ionic current and EOF 


Figure 2a and b show typical variations of the ionic current and flow rate as a function 
of applied voltage, for a negatively-charged pore. Linear fits about Vo = 0 are shown by 
the dashed lines. For positive biases, the magnitude of the current is less than would 
be expected from the linear relationship, while the magnitude of EOF is greater. The 
reverse is true for negative biases. This behaviour we define as rectification, which can 
be quantified by defining rectification ratios as follows: 


n = 

I(-1V) 

n+iv) 

r Q = 

Q(-iv) 

Q(+iv) 


( 6 ) 

(7) 


where the definition has been chosen to be consistent with previous literature. For a 
negatively-charged pore, therefore, rj > 1 , while 0 < cq < 1. 

Figure 2c shows the current and flow rectification ratios as a function of surface 
charge, for different values of taper angle. There are three main observations from these 
results. First, current and flow rectification behave in the opposite sense, as observed 
experimentally m uni- Second, the degree of rectification depends on the surface 
charge, with more highly-charged pores exhibiting greater rectification. A non-zero 
surface charge is therefore a necessary condition for rectification. The final observation 
is that rectification additionally requires an asymmetry such as that introduced by a 
uon-zero taper angle. The greater the pore’s asymmetry, the greater the rectification. 
The aim of this paper is to develop an understanding of these observations. 








(a) 
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Figure 2. Rectification of ionic current and electroosmotic flow, (a) Squares indicate 
ionic current, and the dashed line is a fit about Vq = 0 V. The current grows sublinearly 
at positive voltages, and superlinearly at negative voltages. Here, a = 0.03 rad, 
and Co = 100 mM. (b) A plot of the flow rate for the same pore, showing the 
rectification behaviour in the opposite sense, (c) The rectification ratio, defined as 
ry = |/(— IV)/I(+1V)\ for current (filled squares) and vq = \Q(—\V)/Q(+1V)\ for 
the flow (open diamonds), is plotted as as function of surface charge density for different 
taper angles (measured in radians). Rectification increases with surface charge, and 
also with taper angle. The conditions required for rectification therefore are a non-zero 
surface charge as well as a non-zero taper angle. 
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3.2. Comparison with analytic infinite cylinder model 
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Figure 3. Simulated flux values through the pore, (a-c) The simulated current, 
flow rate, and transference through a cylindrical pore (symbols), and the analytic 
calculations of the same quantities through an infinite cylinder of equal radius (solid 
black lines). Quantities at positive voltages are represented by filled squares, and those 
at negative voltages by open triangles. The transference for both cations (cyan) and 
anions (purple) are shown; the sum of the transference is always 1. Vertical dotted 
lines indicate the limit of validity of the Debye-Hiickel approximation used to derive 
the analytic results, (d f) The same as above, but for a cone of a = 0.03 rad. In this 
case, the magnitudes of the quantites are not preserved under voltage reversal. The 
analytic calculations are now carried out for an infinite cylinder of the same effective 
radius (equal to the average radius of the cone). 


The factors which influence the current and flow can be understood by considering 
a simple model consisting of an infinite cylindrical channel under a constant applied 
electric held. For such a geometry, standard analytic solutions for the ion distribution 
and EOF velocity are available under the Debye-Hiickel approximation \efi/{ksT)\ < 1 
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m- It is therefore possible to analytically determine the integrals in equations [4] and 
[5j For small surface charges, we obtain the following linearized expressions for ionic 
current and flow rate (derivation in Supplementary Information): 


/ = na 2 Den 2 E z 


( 8 ) 



(9) 


where a is the cylinder radius, E z the axial electric field, k = y/2 N,a^ 2 Cq/( k B Te 0 e r ) is 


the inverse Debye length, e = eoe r , and I n is the n th -order modified Bessel function of 
the first kind. The inverse Debye length k oc y/co is a parameter which characterizes 
the salt concentration, and in the rest of this paper the ionic strength will be referred 
to in terms of k. 

In addition, an important quantity to consider is the transference number 7), defined 
as the proportion of total current carried by an ionic species v. T t = I i /I . The linearized 
cation transference number is given by (Supplementary Information): 



( 10 ) 


These equations show that at small surface charges, / is independent of the surface 
charge, while Q and T K + are linear in a s . Figure 3 compares simulated values of /, Q, 
and T k + with these analytic expressions. We observe that the linearized infinite cylinder 
model is an excellent approximation for the finite cylinder at small surface charges, as 


expected (Figure 3a-c). As surface charge is increased, the current does not vary to first 


order, while the flow rate increases linearly. The cationic transference also increases, 
showing clearly the increasing permselectivity of the pore. 

Once a taper angle is introduced, the infinite cylinder model still captures the 
qualitative trends in /, Q, and T K + 1 with one important exception: the simulated 
results show that the magnitude of these quantities is no longer the same upon voltage 
reversal (Figure 3d-f). At negative voltages, the magnitude of / is increased compared 
to at positive voltages, while Q and T^+ are smaller. 
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Equations ii and [lO] suggest that the only quantities that can change upon 
voltage reversal are electric field and salt concentration, since all other quantities 
remain constant. Because the equations do not determine E z and k, requiring instead 
that they are specified by other means, they cannot explicitly capture the rectification 
behaviour. However, the main conclusion that can be drawn is the fact that the change 
in magnitudes of /, Q, and T K + upon voltage reversal are due to a change in electric 
field and salt concentration. 


3.3. Relationship between electric field and salt concentration 

In order to investigate the relationship between electric field and salt concentration, we 
begin by calculating the average value of these quantities over the entire volume of the 
pore. Figure 4a and b show the average electric field magnitude and salt concentration 
(measured in terms of the inverse Debye length k) inside the conical pore, as a function 
of surface charge, for voltages of Vo = ±1 V. Here, the difference in behaviour due 
to the sign of the applied voltage becomes apparent. In Figure 4a, we see that for 
a given surface charge, positive voltages result in high electric fields, while negative 
voltages result in low fields. The opposite effect is true for the salt concentration: 
positive voltages result in a lower salt concentration compared to negative voltages. 
The variations in electric field and salt concentration are thus anti-correlated. 

We can investigate the dependence of electric field on concentration in more detail 
by considering the relationship between the magnitude of flow rate and transference, 
as plotted in Figure 4c. This shows that the flow is a linear function of transference. 
The intuition behind this observation is that, for symmetric electrolytes with equal 
mobilities, cationic transference is essentially a measure of the prevalence of cations 
over anions, and hence the degree of momentum imbalance in the fluid. Without an 
imbalance, i.e. for T K + = Tqi- = 0.5, there can be no net force and hence no EOF. By 



Electric field (10 6 V/m) 


12 


(a) 




cone 


cone 



Figure 4. Relationship between electric field and salt concentration, (a) The average 
electric field magnitude inside a conical pore with a = 0.03 rad as a function of surface 
charge density. For positive biases, as charge is increased, the field magnitude increases, 
while at negative voltages, the field decreases in magnitude, (b) The average salt 
concentration (measured in terms of k, the inverse Debye length) exhibits the opposite 
behaviour. As charge is increased, concentration decreases at positive voltages, but 
increases at negative voltages, (c) The magnitude of the flow rate is a linear function 
of transference, and the slopes of the relationship at positive and negative voltages 
are roughly equal. This sets a constraint for the relationship between the electric field 
and salt concentration, (d) The electric field magnitude and concentration inside the 
pore are anti-correlated. Negative voltages tend to drive the pore into a low held, high 
concentration state, while positive voltages result in a high held, low concentration 
state. The solid line is the analytic relationship expressed in equation [j~2| 
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eliminating a s between equations [9] and 10, one obtains 

1 2 {na 


Q = 


7ia 2 e 


I ( k B T \ (na Io(na) _ i A IAkCL) 
iT ' \ D J 2 /,(m) 1 / 


na)E z 


T - - 2 


( 11 ) 


ekgT 1 V D ) \ 2 Ii(kcl) 

The observed linear relationship in Figure 4c, and the nearly equal slopes for both 


positive and negative applied voltages, indicate that the prefactor in equation 11 


is 


approximately constant. It is important to point out that a priori , there is no reason to 
expect that the slopes should be equal upon voltage reversal. The consequence of this 
observation is that the variation in electric field and salt concentration is constrained to 
obey 


7ra 2 e 


hina) 


e~ht f k B T \ ( na Iojua) _ -i A 

ek B T ' V D ) 2 /iH 


( na)E z ~ C, 


( 12 ) 


where C = Q/T K +. Figure 4d shows the electric field against salt concentration, and 


the analytic relationship given by equation 12 This plot clearly shows the constraint 
on the variation of electric field and salt concentration inside the pore: low fields are 
correlated with high concentrations, and vice versa. 

An intuitive understanding of the anti-correlated behaviour of electric field and salt 
concentration can be obtained by considering Ohm’s law, which states that J = oE , 
where cr is the conductivity of the solution. For a constant current, as the conductivity 
reduces, the electric field increases. This would result in a constraint E oc /t -2 , i.e. the 
same as that expressed by equation [8] However, it is important to note that under 


rectification conditions, the current is no longer constant. Since equation [12] takes into 
account varying current as well as EOF, it is thus a more sophisticated constraint than 
Ohm’s law. 


3-4- The relationship between ionic current and EOF 


Under situations where the constraint given by equation 12 is valid, the ionic current 


and flow rate can be written down purely in terms of the salt concentration. 

kIAko) 


I = AC 


h (ko) 


^ f na Inina) 

1 + S- ——- - 1 

1 2 h{na) 


(13) 
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Figure 5. The variation of current and EOF with salt concentration, (a) The current 
is an increasing function of salt concentration, as described by equation |13| Values 
plotted here are in arbitrary units, (b) The flow rate is a decreasing function of salt 
concentration, as described by equation [l4| 

Q = Bc\ 

Av 

where the constants A = Defi/(ksTa), B = —<r s e/(e/csTa), and S = e(ksT) 2 / (De 2 fi). 
Only one (E, k) pair is needed to set the constant C, and hence in these expressions the 
functional dependance of the electric field is removed completely. 

Figure 5 shows the form of these two equations. Whereas the ionic current is 
an increasing function of salt concentration, the flow rate decreases as concentration 
increases. A change in salt concentration thus leads to variation of current and flow in 
opposite directions. 

Whilst we have described the variations of electric field and concentration in detail, 
we have made no mention yet of the cause of these variations: in particular, what is 
required is a mechanism which results in different concentrations inside the pore as a 
function of voltage. 


(14) 
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Figure 6. Concentration polarization in the nanopore, (a) The variation of salt 
concentration along the axis of the pore, for voltages of ±1 V applied from left to 
right; vertical lines indicate the extent of the pore. An ionic current flowing through 
a permselective pore results in a depletion region upstream, and an enrichment region 
downstream of the pore. For a cylinder, the polarization is symmetric about the 
centre of the pore. Different lines correspond to different values of surface charge 
density, which is increased in equal intervals from 0 to -40 mC/nr 2 . (b) When the pore 
contains a taper angle of 0.03 rad, the CP is asymmetric with respect to the centre 
of the pore. At positive voltages, the upstream depletion region is extended inside 
the pore compared to the enrichment region outside, (c-d) The spatial variation in 
concentration near the same conical pore (with a s = —0.03 C/m 2 ), showing clearly 
the depletion at positive voltages and enrichment at negative voltages. 
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3.5. Concentration polarization and symmetry breaking 

As discussed earlier, an ionic current flowing through a permselective pore results in 
a depletion region upstream, and an enrichment region downstream of the pore, an 
effect called concentration polarization (CP). CP arises due to the requirements of flux 
balance. Within the permselective region, the majority of the cationic flux is due to 
electromigration; however, upstream of this region, the same flux is distributed between 
a smaller electromigration component, and a diffusive component towards the pore. 
The diffusive flux results from a concentration gradient which decreases approaching 
the pore, which characterizes the upstream depletion region. 

Conversely, downstream of the pore the flux is distributed between an 
electromigration component and a diffusive component away from the pore, which results 
from a concentration gradient decreasing away from the pore. This characterizes the 
downstream enrichment region. 

Electroneutrality requires that the anionic concentrations are depleted and enriched 
in the same way as the cations, and thus the variation applies to the overall salt 
concentration, not just that of a single species. 

In a cylindrical pore, this polarization is symmetric with respect to the centre of the 
pore (Figure 6a). However, by introducing a taper angle, the polarization can be made 
asymmetric with respect to the direction of current flow: depletion and enrichment 
regions inside the pore become much larger than the corresponding regions outside 
(Figure 6b). Hence, the entire pore can be driven into a depleted, low-concentration 
state under positive voltages, and a high-concentration state under negative voltages 
(Figures 6c and d). This asymmetric behaviour is the origin of the rectification effects. 

The final picture we have for current and flow rectification can therefore be 
summarized as follows: in a permselective, conical pore, a positive voltage leads to CP 
which reduces the salt concentration inside the pore. This reduction in salt concentration 
puts the pore in a low-conductivity state, and hence a low ionic current flows. At 
the same time, the low concentration is correlated with high electric fields and a high 
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cationic transference number, as cations form a larger proportion of the total ionic 
concentration. The high electric fields and prevalence of cations over anions result in 
large electroosmotic flows. The converse, when a negative voltage is applied, results in 
a high conductivity, low transference, and low flow state for the pore. 

Because CP is absent in a non-permselective pore, the presence of surface charge is a 
necessary condition for rectification. In order to obtain asymmetric CP, an asymmetric 
pore is required, which is the second necessary condition. Our model thus correctly 
predicts the conditions required for rectification as presented in Figure^. 

4. Discussion 

We have presented a unified description of ionic current and electroosmotic flow 
rectification. The study has been guided by analytic relations derived from an infinite 
cylinder model, and is therefore valid for small surface charges, small taper angles, and 
relatively large aspect ratio pores (which are long compared to their width). These 
requirements have motivated the conditions we have used: surface charges of 0 to -0.03 
C/m 2 compare well with literature values of charges on glass and silica surfaces [5j, and 
are small enough that the Debye-Hiickel approximation is valid for a reasonable range 
of values. We have presented results for taper angles between 0-0.03 radians, which 
corresponds to a maximum total opening angle of ~ 3°. Experimentally, conical glass 
nanopores can have opening angles up to ~ 10° [IS]. Finally, our model describes conical 
pores with a length of around 100 nm, and can be best applied to systems such as conical 
glass nanopores D2J- Since typical solid-state nanopores are made in membranes with 
thicknesses on the order of a few tens of nm [26], we expect effects of the membrane 
outside these pores to play a more important role mm, and hence our description of 
the pore as a long cylinder will require modification when applied to these systems. 

As the parameters vary outside these ranges, we expect a departure from linearity. 
This is already exhibited in the increasing current, and sublinear growth in flow and 
transference (e.g. Figure 3a-c) at high surface charges. This behaviour is due to the 
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salt concentration inside the pore being effectively dominated by the ions inside the 
double layer. As the pore becomes more highly charged, electroneutrality requries that 
the total charge in the double layer balances the total charge on the walls, and so the 
salt concentration inside the pore slowly increases beyond even the bulk values. This 
leads to higher ionic current and a slower growth in the flow, as expected. 

Although we have presented results for fixed pore radius and salt concentration, 
an increase in either radius or salt concentration will diminish rectification effects, 
as observed experimentally (e.g. TO- This can be attributed to the decreased 
permselectivity of the pore in both cases. 

In the Supplementary Information, we show results for charge densities up to -0.08 
C/m 2 and taper angles of a = 0.5 rad, as well as the effects of varying the pore length 
and radius. In this case, although the linear approximations break down, the results 
preserve their qualitative properties, and we believe that as long as the system remains 
free from instabilities, such as those associated with overlimiting current behaviour m, 
our conclusions will remain valid. 

As mentioned earlier, it is intuitive that flow rate and transference should be 
correlated, as the cationic transference measures the charge imbalance which leads to 
a net force in the first place. However, the observation that the relationship between 
the two quantities is linear (Figure 4c) is a stronger statement, which allows us to write 
down the constraint given by equation [l2j The origin and generality of this constraint 
will be investigated in future work. 

Finally, although in our studies the asymmetry which is responsible for rectification 
is geometric, we expect that any other property which breaks the (—z —> z) symmetry 
of the system will lead to similar EOF rectification effects. An example would be a 
spatially-varying surface charge, which has already been demonstrated to lead to current 
rectification in cylindrical pores |28| . as well as the closely-related effect of osmotic flow 
rectification [29] . 
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5. Conclusions 

We have carried out an investigation into the rectification of electroosmotic flows in a 
conical nanopore. Like current rectification, this effect requires the pore to possess both 
permselectivity and geometrical asymmetry. The flux of current and water through the 
pore are determined by the electric field and salt concentration inside the pore, two 
quantities which we have found to be anti-correlated. Permselectivity and asymmetry 
result in the pore exhibiting concentration polarization: when a positive voltage is 
applied at the wide end of the pore, the pore is driven into a low concentration state 
associated with low current and high electric fields, transference, and EOF. A negative 
voltage drives the pore into a high concentration state, which results in high currents 
and low electric fields, transference, and EOF. The rectification of current and flow is 
therefore intricately linked. 
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